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We sense the motion of a trapped atomic ion using a sequence of state-dependent ultrafast mo¬ 
mentum kicks. We use this atom interferometer to characterize a nearly-pure quantum state with 
n — 1 phonon and accurately measure thermal states ranging from near the zero-point energy to 
h ~ 10 4 , with the possibility of extending at least 100 times higher in energy. The complete energy 
range of this method spans from the ground state to far outside of the Lamb-Dicke regime, where 
atomic motion is greater than the optical wavelength. Apart from thermometry, these interfero¬ 
metric techniques are useful for characterizing ultrafast entangling gates between multiple trapped 
ions. 


There is great interest in the ultrafast quantum control 
of trapped ions, including the preparation of nonclassi- 
cal states of motion |Tj, entangling quantum logic gates 
0 0 , and ground state cooling on time-scales shorter 
than the period of ion harmonic motion [4] . Experi¬ 
mental demonstrations of ultrafast trapped ion control 
include single qubit rotations m Ei and spin-motion en¬ 
tanglement 00 . These operations can be many orders 
of magnitude faster than conventional techniques, less 
sensitive to noise, and scalable to large numbers of ions 
where the couplings occur through local modes of motion 
0. Moreover, these ultrafast techniques work far outside 
the Lamb-Dicke regime, where the extent of the atomic 
motion is greater than the optical wavelength 0 . 

Ultrafast sensing of atomic motion allows measure¬ 
ments over a wide range of energies, from the zero-point 
(average phonon ocupation number n = 0) to poten¬ 
tially above room-temperature (n ~ 10 6 for typical ion 
traps). Ultrafast partial state tomography (defined later 
in this letter) on thermal states improves upon the dy¬ 
namic range achieved with thermometry using dark reso¬ 
nances m- It also complements conventional methods of 
thermometry, including measurements of the motionally- 
induced upper and lower sideband asymmetries m and 
the thermal suppression of induced transitions (Debye- 
Waller factors) [T2j. However, both of these other meth¬ 
ods break down when the atomic motion is outside of the 
Lamb-Dicke regime, typically around n > 10. Measuring 
the entire Doppler-broadened envelope of sidebands [13] 
provides a more general measurement of thermal states, 
but can be difficult due to the bandwidth required to 
excite multiple sidebands. Here we use ultrafast tech¬ 
niques for accurate thermometry of ion motion ranging 
from n ~ 0.1 to n ~ 10 4 and show how this method 
extends to higher energies. We also measure particular 
quantum states through more complete motional tomog¬ 
raphy. 

In this experiment, we trap a m Yb + ion in a linear 
radio frequency Paul trap and probe the motion along 
a single radial mode of motion with secular trap fre¬ 
quency wt/ 27t « 1 MHz, as described in Ref. 0. The 
|.F = 0, to/ = 0) = ||) and |F = 1 ,to/ = 0) = ||) hyper- 


fine levels of the 2 Si/ 2 electronic ground state are used 
as the qubit, or effective spin, and are separated by the 
splitting = 12.642815 GHz. The ion is laser-cooled 

to near the Doppler limit using the 2 Si/ 2 to 2 Pi/ 2 transi¬ 
tion at a wavelength of 369.5 nm, and optically pumped 
via the 2 Pi/ 2 , F = l levels to the ||) level during initial 
state preparation. Qubit state detection is performed by 
collecting state-dependent fluorescence [14]. The qubit 
state in these experiments is detected with an efficiency 
above 0.997 using an imaging objective with 0.6 numeri¬ 
cal aperture and a photomultiplier tube m- 

We create a spin-dependent kick (SDK) by shaping in¬ 
dividual pulses extracted from a mode-locked laser with 
center wavelength 2n/k ss 355 nm, pulse duration r ~ 10 
ps, and repetition rate of f rep = 118 MHz. A series of 
optical delay paths in the form of three sequential Mach- 
Zehnder interferometers shape the pulse, dividing it into 
eight sub-pulses that are spaced to collectively flip the 
qubit spin while generating spin-dependence in the mo¬ 
mentum transfer 0 . In addition to setting proper de¬ 
lay lengths, we adjust the pulse energy to give a com¬ 
plete spin flip 0 . The pulse is applied to the ion 
in a counter-propagating geometry, ideally creating the 
evolution operator Usdk = 'D[ir)\a+ + T>[—ir]\a- TB] , 
where a± are the qubit raising and lowering operators. 
The displacement operator T>[±irj\ imparts momentum 
A p = ±2hk = ±r/po, where po = y/2MhuJt is the zero- 
point spread of momentum in the trap, M is the atomic 
mass, and p = 0.2 is the Lamb-Dicke parameter associ¬ 
ated with this momentum transfer. The impulsive SDK 
operation occurs on a time scale much faster than the 
trap period (r <C 1/w*), and the spin population transfer 
from ||) to |t) is measured to have a fidelity of 0.993(2) 
0 . Because each SDK operation provides a momentum 
kick and flips the spin, immediately applying a second 
SDK would simply undo the first. However, by waiting 
one half of the trap period between SDKs, we can con¬ 
catenate N individual kicks to create a larger effective 
SDK with A p = ±2 Nhk = zLNr/po. 

Techniques dealing with Ramsey spectroscopy on 
states coherently displaced by spin dependent forces have 
been demonstrated in creating Schrodinger cat states m 
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and measuring spin dephasing in 2D ion crystals [18] . In 
this experiment, we create an interferometer to sense mo¬ 
tion by applying two sets of N SDK operations within a 
Ramsey experiment on the qubit levels with time dura¬ 
tion T. First the ion is prepared in a coherent superposi¬ 
tion of 14,) and |t) by applying a near-resonant microwave 
7 t/ 2 pulse of duration r M . Following the first set of N 
SDKs, the ion evolves for a time 9/u> t and then a second 
set of N SDKs is applied. Finally, another ir/2 pulse with 
the same duration and tuning drives the qubit to close the 
Ramsey interferometer. This sequence is diagrammed in 
Fig. 0c By scanning the microwave detuning S <C 1 /r M 
from resonance, we observe sets of Ramsey fringes with 
phase <fi = 6T that chronicle the ion motion (shown in 

Fig. 0> and So- 

For a pure initial state |^“) i = |4-) |a), where a is a 
coherent state of the ion motion, the state following the 
Ramsey experiment is 

|*“> = ^ 1 {\i)+ie- i< ‘ > \\))\{a + iNr 1 )e- ie -iNr]) + 

*e“* 7 (|t) + |4-)) |(a - iNr])e~ l6 + iNrj)], 

( 1 ) 


where 7 = iV 77 [Re(a)(l — cos 9) — Im(a) sin 9\. 

Given an arbitrary initial state of motion in phase 
space described by the Glauber P-distribution pa Ea, 
the final density matrix is p = f P(a) |H/ a ) (\I/“| d 2 a. The 
probability of measuring the state spin-up after the Ram¬ 
sey experiment is therefore 

S{9,N;4i) = (f|p|t) 

= ^ ^ J P{a)e~ i( ' Nr, ' > i 1 ” 0080 ) cos(47 — cj))d 2 a. (2) 

Two types of motional state that are readily accessi¬ 
ble in the laboratory are thermal states and small Fock 
states. First we discuss ultrafast partial state tomogra¬ 
phy to determine the average phonon number in a ther¬ 
mal state. Then we extend this method to create a nearly 
complete map of the motion of an n = 1 Fock state in 
phase space, showing clear nonclassical signatures. 

For an ion prepared in a thermal state with 
mean phonon number n and P-function Ptherm{oc ) = 
_L e -M 2 /™, Eq. [ 2 ] yields an expected Ramsey fringe pat¬ 
tern 

S th erm(9,N 1 + ie-^^+W 1 - 0 " 9 W- (3) 

The fringe contrast has periodic revival peaks at 9 = 
27tto, where to is a positive integer. For a hot ion where 
n l/(Nr)) 2 , these revivals in contrast become narrow 
and approximately Gaussian with full width at half max¬ 
imum FWHM= 0.83 /(Nrj\/n). With N = 1, we measure 
the Ramsey fringe contrast as a function of 9 for a variety 
of initial thermal states of motion, and fit the contrast 
revival peaks to Eq. 3 to determine the average phonon 
number n of the thermal state [2T1 - |23| . In the fit, we al¬ 
low the peak Ramsey contrast at 9 = 2i rm to be less than 
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FIG. 1. (a) Timeline of a single experiment, where a full 

SDK set is made of N single SDKs, (b) Phase space diagram 
of an initial state (|4-) + |T)) |«) evolving under two sets of 
SDKs separated by time delay 9/uJt, where |a) is a coherent 
state of motion, (c) Typical Ramsey fringes as a function of 
microwave frequency detuning S. These two plots correspond 
to the points 9 = 0 and 6 = 7r/2 of an initial thermal state 
(N = 1 for the data shown). 


unity in order to parametrize imperfect fidelity of the 
SDK operations. This reduction in fidelity is mainly at¬ 
tributed to variations in the Raman beam intensity over 
the spatial extent of the ion wave packet, and becomes 
apparent at high n. This does not significantly affect the 
width of the contrast revival peak or the accuracy of the 
thermometer, and can be mended by widening the beam 
waist. 

Ramsey contrast revival lineshapes are measured in ex¬ 
periments spanning over five orders of magnitude in n, 
as shown in Figs. 2a-c. Figure 01 shows these mea¬ 
surements plotted versus the expected value of n from 
theory and other measurements. The figure is broken 
into three regions according to the manner in which the 
motional state is prepared and calibrated before mea¬ 
surement of the contrast revival lineshapes. Low energy 
thermal states (h < 10) are generated by first sideband¬ 
cooling the ion to its zero point motion and then allowing 
the ion to weakly heat in the trap by known amounts. In 
this regime, we compare ultrafast interferometric mea¬ 
surements of n (shown in Fig. 0) to values extracted 
from measured sideband asymmetries [241 . The deviation 
of the two measurements are shown in the red section of 
Fig. 0. 

For thermal states 10 < n 4$ 150, the ion is prepared by 
Doppler cooling with various frequency detunings from 
resonance. Ultrafast measurements in this regime are 
shown in Fig. 0. Each of these measurements is com¬ 
pared to the predicted value of n from Doppler cooling 
theory [25] and plotted against each other in the green 
section of Fig. 0. As a check on the expected values of n 
in this range, we also measure the Debye-Waller suppres- 





















3 




9/2tt 



0/2n 



FIG. 2. Ultrafast sensing measurements of n (with 
N = 1): (a) Sampling of Ramsey revival contrast lineshapes 
with initial states prepared by resolved sideband cooling to 
the ground state and subsequent heating. Data is fit to 
Stherm(d, 1; </>)• The amplitude of each fit is a free parameter 
to account for SDK infidelity (also done in b and c). This 
does not significantly affect the width of the peak, which is 
used to determine h. (b) Sampling of Ramsey revival con¬ 
trast lineshapes with initial states prepared by Doppler cool¬ 
ing only, with n varied by changing the cooling beam de¬ 
tuning. (c) Sampling of Ramsey revival contrast lineshapes 
with initial states prepared by inducing a high heating rate 
with white noise applied to a trap electrode, (d) Measure¬ 
ments of n versus predicted values. There are three regimes 
of thermal state preparation-red being sideband-cooling then 
heating (see (a)), green being Doppler cooling with different 
detunings (see (b)) and orange being heating with applied 
noise (see (c)). 


sion of Rabi flopping transitions between the ion qubit 
states [24] for several cases, resulting in expected values 
consistent with Doppler theory. 

Hot thermal states are prepared by inducing a high 
heating rate with a noisy electrical potential to a trap 
electrode for varied amounts of time after Doppler cool¬ 
ing. The ultrafast measurements of these states are 
shown in Fig. w Measurements in this regime are com¬ 
pared to a predicted n given by the equation n = 

HU, where e is the ion charge, and Sv(wt) (V 2 /Hz) is 
the applied power noise spectral density of the electric- 
potential, which is white over the measurement band¬ 
width. The effective distance d of the electrode to the 
ion is calibrated by applying a static potential offset to 
the same electrode and observing the resulting displace¬ 
ment of the ion in space. The predicted and measured 
values for this regime are plotted against each other in 
the orange region of Fig. § 1 - 

We next perform more complete tomography of a 
nearly pure quantum state of motion by extracting the 
characteristic function 

Xw(o 0 = e - |Q|2/2 J P(/3)e Mm(a/ 3 * ) d 2 / 3. (4) 

This quasiprobability distribution contains all the infor¬ 
mation about the quantum state and is the Fourier trans¬ 
form of the better-known Wigner distribution [251 [23- 

In terms of the observable S(0,N;</>), xw(a) is given 
by 

Re[ Xw (a)]=2S(0,iV;O)-l (5) 

Im[ Xw (a)]=2S(0,JV;|)-l, 

where a = 2Nrj[sin9 + i( 1 — cos#)]. Scanning 6 and N 
while measuring S(9, IV; </>) maps the characteristic func¬ 
tion over rings in phase space, shown in Fig. §i. In 
order to scan the negative imaginary part of a, we can 
change the direction of the initial momentum kick asso¬ 
ciated with the spin flip operators by shifting the relative 
optical phase of the counter-propagating beams by ir [7]. 
These reversed kicks can be thought of as effectively flip¬ 
ping the sign of 77 , and for simplicity, we represent them 
here by negative values of N. 

We measure the characteristic function Xw( a ) of the 
ion in the n=l Fock state, prepared by sideband cooling 
to the ground state and transferring population to the 
n=l state through application of a blue sideband oper¬ 
ation 9j. To have a grid that spans the domain of the 
state, we scan around 16 rings in phase space set by ±IV, 
where N = 1, 2,3,4, 5, 6 , 8 ,10. Two of the 16 rings along 
which we measure are highlighted in Fig. §b and plots 
of S(9,N; 0) versus 9 along those two rings are shown 
in Fig. §>. Notice in Fig. [3]o that the larger SDK set 
(N = 5) separates the interferometer enough to see the 
oscillation of the motional distribution, while the smaller 
SDK set does not. Mapping along all 16 curves gives a 
nearly complete motional state map. The real part of 
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FIG. 3. (a) Points in phase space accessible in our tomographic measurements. The radius of each circle (2Nr]) is set by the 
number of kicks N, and the angular position on each circle is set by the SDK delay 6. The sign of N represents the direction of 
the initial momentum kicks associated with the spin flip operators, (b) Sample of measurements of the Ramsey fringe at if> = 0 
for a nominal n — 1 Fock state, using two sets of kicks with N — 5 (red) and N = —2 (blue) and scanning the delay 9. The 
coordinates of these particular scans in phase space are highlighed in a. (c) Motional state tomography of an ion prepared in 
the n = 1 Fock state. In ascending order: the value-colored data points of Re[\vv(aO] taken on 16 rings in phase space set by 
±N where N = 1,2, 3, 4, 5,6, 8,10, the interpolated data mapped in contour, and a 3D interpolation of the data, (d) Theory 
prediction for a Fock state with n = 1. 


the characteristic function is shown in Fig. H alongside 
the corresponding model of Re[xw(oO] for an n = 1 Fock 
state in Fig. W- The negative values of the character¬ 
istic quasiprobability function highlight the nonclassical 
nature of the motional state of the ion. 

These ultrafast tomographic techniques are capable of 
measuring motional energies far beyond the data pre¬ 
sented here, which was limited to n ~ 10 1 2 3 4 5 because of 
re-cooling issues during state preparation. In the ex¬ 
periment, we scan the interferometric angular delay 6 
in steps set by the repetition rate of the laser, giving a 
resolution of u>t/ frep ~ 50 rnrad. For lineshapes narrower 
than this laser repetition rate limit, we scan 6 by chang¬ 
ing the trap frequency u>t though accurate control of the 
trap rf drive voltage. With fine drive voltage control, we 
can achieve a resolution in 9 of 0.1 mrad, which would 
correspond to a Ramsey revival linewidth from a ther¬ 
mal state with h ~ 10 9 . Other factors also come into 
play when measuring such high-energy states: First, the 
spatial extent of motion swells beyond the laser beam 
waist. At fi = 10 6 for instance, or equivalent tempera¬ 


ture T = hu}t.n/kB = 80K, the ion would experience a 
significant gradient in the Rabi frequency across a beam 
with a 3/rm waist. A second factor is the decreased de¬ 
tection fluorescence due to larger Doppler shifts at these 
energies. The detection fluorescence at n = 10 6 would be 
reduced by a factor of ~ 10 3 from a cold ion [25|. Finally, 
when measuring these very narrow lineshapes, instabili¬ 
ties in the trap frequency LOt and laser repetition rate f rep 
would have to be sufficiently stable over the measurement 
time. At n = 10 6 , this would require a fractional stability 
from both the trap frequency and laser repetition rate of 
better than 0.1%. These factors put ultrafast interfero¬ 
metric measurements of n > 10 6 neither fundamentally 
nor technically beyond reach. 
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